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Abstcast 

A loosely coupled two-phase vacuum water plume 
model has been developed. This model consists of a 
source flow model to describe the expansion of water 
vapor, and the Lagrangian equations of motion for 
particle trajectories. Gas/Particle interaction is modeled 
through the drag force induced by the relative velocities. 
Particles are assumed traveling along streamlines. The 
equations of motion are integrated to obtain particle 
velocity along the streamline. This model has been 
used to predict the mass flux in a 5 meter radius 
hemispherical domain resulting from the burst of a 
water jet of l.S mm in diameter, mass flow rate of 24.2 
g/s, and stagnation pressure of 21.0 psia, which is the 
nominal Orbiter water dump condition. The result is 
compared with an empirical water plume model, deduced 
from a video image of the STS-29 water dump. To 
further improve the model, work has begun to 
numerically simulate the bubble formation and bursting 
present in a liquid stream injected into a vacuum. The 
technique of smoothed particle hydrodynamics was used 
to formulate this simulation. A status and results of 
die on-going effort are presented and compared to results 
from the literature. 


Introduction 

The Orbiter dumps a combination of waste and 
condensate water from a port-side, heated, knife-edge 
nozzle, approximately every 3 days. The narrow jet of 
water injected into vacuum quickly bursts into a 
disperse plume of vapor, water droplets, and ice 
particles. The expanding plume may make contact with 
surfaces, and therefore stick and leave solid contents as 
permanent deposits. The optical properties of the 
impacted surfaces may be affected to a degree that 
inhibits their proper function. Phenomena that occur 
when a liquid stream is injected into a vacuum 
environment have been studied by Fuchs and Legge 1 and 
by Muntz and Dixon 2 . In general, as a liquid stream is 
discharged into a vacuum, due to the sudden drop of 
pressure, the stream becomes superheated and vapor 
bubbles may form inside the stream. As the bubble 
continues to grow, it eventually reaches a critical size 
which causes the jet to burst into water droplets and ice 


fragments. The fragments form a conically shaped 
plume. The burst distance and the characteristic cone 
angle are affected by the initial stream velocity, 
temperature and diameter. The effect of these parameters 
on the burst distance and the characteristic cone angle 
were studied. The results of these parametric studies are 
presented in the following section. For the purpose of 
predicting the extent of contamination due to the Orbiter 
water dump, an engineering model has been developed 
to model the fast expanding gas vapor and the conical 
cloud of ice fragments. This model consists of a source 
flow model to describe the vapor phase of the plume, 
and the Lagrangian equations of motion for particle 
trajectories. The interaction between the particle and 
gas is modeled through only the drag induced by the 
relative velocities 3 . Particle number density along the 
plume centerline was deduced from a vacuum venting 
test conducted at Arnold Engineering Development 
Center (AEDC) 4 - 5 . A normal distribution function is 
used to describe the particle number density in the polar 
direction. The variance of this distribution may vary and 
two values of distribution variance based on the plume 
cone half angle have been examined. Given a location 
with respect to the nozzle exit, this model can predict 
mass flux at the specified location due to the dump. 
This paper will provide a detailed description of the 
model and the numerical procedure used. This model has 
been used to generate the mass flux in a 5 meter radius 
hemispherical domain resulting from the burst of a 
water jet of 1.5 mm in diameter, mass flow rate of 24.2 
g/s and stagnation pressure of the jet is equal to 21.0 
psia, corresponding to the nominal Orbiter water dump 
condition. The result is compared with the result 
obtained from a purely empirical model, deduced from a 
video image of the STS-29 water dump 6 . Model 
improvements via a new numerical simulation will then 
be discussed. 


Cavitation of Liquid Streams in a Vacuum 


The equation of motion for the radius R b of a 
bubble in a viscous liquid as a function of time is 


R b R b +\ R b 


AP 2a 4/jR b 
Pi Pl R b Pl R b 


( 1 ) 
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where a is the surface tension, p,is the density of the 
liquid stream, and p is the viscosity 7 . AP is the 
pressure difference driving the bubble expansion and can 
be expressed as 


where p vb is the vapor pressure of the liquid at the 
bubble wall temperature. The vapor pressure is given 
by the Clausius-Clapeyron equation 

Pv =PR*xp(-^r) (3) 

kT s 

where p R is a constant obtained empirically, m, is the 
mass of a liquid molecule, l y is the heat of vaporization, 
k is the Boltzmann constant, and F, is the surface 
temperature. The reduction of the surface temperature 
with time of a water stream due to evaporation as given 
in Fuchs and Legge 1 is 


T 0 ~T s (t) = 



( 4 ) 


where a is the thermal conductivity, T 0 is the initial 
stream temperature, and T g (t) is the stream surface 
temperature at time The time-dependent surface heat 
transfer for vaporization is 


p b = p s + a/R 0 + 2a/R b 

where p, is surface pressure of the jet stream, R 0 is the 
radius of the stream, and R b is the radius of the bubble. 
If we assume there are no temperature gradients in the 
stream, the internal bubble pressure must be the vapor 
pressure and p, = p v /2, or p s = 0 few a very cold 
surface. The bubble must have at least a radius of R b as 
solved for in equation (7), if the bubble is to begin 
growing. The equilibrium bubble radius, or minimum 
radius that will permit growth is 

Rbe = (p v /2-a/Ro) (8) 

or, for a odd stream surface. 




2a 

(Pv-C/Ro) 


( 9 ) 


It is assumed that bursting takes place when a 
bubble grows to be the size of the stream, and that the 
perpendicular particle velocity is equal to the growth 

rate of the bubble, R b . Therefore, the cone half angle 
of the region containing majority of ice fragments is 
given by 


0= arctemf R^ Vq) (10) 


q(t)srB*XL, |i*El 
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( 5 ) 


where n vs is the vapor number density at temperature 
T s . The definite integral of equation (4) has the 
interesting property that it is dominated by the value of 
integrand near t—t.. A useful approximation, as given 
in Muntz and Dixon 2 , which permits a simple explicit 
expression to be derived, is to assume that q is 
constant at its value at the end of die time interval. It 
follows that the surface temperature can then be 
obtained by solving the following non- linear equation 

To -T,mS^m ( 6 ) 

a 

where c p is the specific heat at a constant pressure. 

Equation (1) is solved by using a third order Runge- 
Kutta integrator 8 . At each time step, equation (6) is 
solved to obtain stream temperature. From Muntz and 
Dixon 2 , in order for a vapor bubble to exist or grow, its 
internal pressure must be at least 


where V 0 is die initial stream velocity. Equation (1) 
has been solved for a 1.5 mm stream at various 
temperatures and velocities. Figure 1 shows the effect 
of stream temperature to the bubble growth rate. The 
result is plotted as Rf/R 0 versus a non-dimensional 
time (tV 0 /R 0 ). Note that the stream is assumed to 
have burst at Rf/R 0 = 1. From Figure 1, we can see 
that the burst distance decreases with increasing stream 
temperature. Figure 2 shows the effect of stream 
velocity cm the burst distance. Again, the result is 
plotted as Rj/R 0 versus the non-dimensional time 
From this figure, we can see that the burst distance 
increases with increasing stream velocity. Figure 3 is a 
composite plot of burst distance versus stream 
temperature at various stream velocities. Again, the 
decrease of burst distance with increasing temperature is 
clearly shown. Figure 4 is a composite plot of burst 
distance versus stream velocity at various stream 
temperatures. Figure 5 is a composite plot of burst 
angle versus temperature at two different velocities. 
From this figure we can see that the burst cone angle 
increases with increasing temperature, and also the cone 
angle decreases with increasing velocity. 
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Plume Flow 


Particle Phase 


After bursting, the liquid stream becomes a cloud of 
gas, liquid droplets and ice fragments. In this work to 
date, we assume all liquid droplets are solidified. 
Therefore, only the gas phase and ice particles are being 
considered in this paper. However, with the SPH 
technique and the SPHINX model, our ul tima te goal is 
the modeling of gas, liquid drops, and solid particles. 

Gas Phase 


The fast expansion of the gas phase is characterized 
by high Mach numbers and velocity; and, therefore, 
almost straight streamlines, which seem to originate at 
the burst point. This type of flow can be describe by a 
source flow model 9 . In such a flow the density at a 
distance r from the burst point is given by 


P_ 

Po 


= ^-(cos 


*JL\Wr-n 

2 9 Um 


( 11 ) 


The particulate model is taken from Holcomb 3 , 
which implements a Lagrangian treatment for the 
particles. The particles are assumed to be solid (no 
mass transfer or size change) and dilute (negligible 
volume in comparison to the gas), but the gas/particle 
system is assumed to be in non-equilibrium, in that 
there are relative velocities between ice particles and the 
expanding gas. The Lagrangian treatment of the 
particles, also referred to as the particle trajectory model, 
utilizes the natural equations of motion of the individual 
particles. 

Particle/gas momentum interaction term accounts 
for the acceleration which the gaseous phase imparts on 
the solid particles. The Lagrangian equations of motion 
for the particles are 


dR p 

dt 


= u r 


(16) 


where p is the stagnation density, 8 is the polar angle, 
8 lim is the limiting angle of the expanding gas, and yis 
the ratio of specific beats. Apis a constant which is 
calculated by mass flow considerations 

A p = --3 /[ 2uUm) (12) 

l‘ im (cos—— — ) 2/(r ~ ]) sinddd 

Jo 2 6 Um 

where ii* is the velocity at sonic condition. u lim is the 
limiting velocity and is given by 

u Um =i^-R T o (13) 

V y- 1 


where R is the specific gas constant and T 0 is the 
stagnation temperature. The velocity at a distance r from 
the burst point is computed from the isentropic 
expansion relation 


U 2 [ 


1_ 1 , dp 

M T (y^T) + 2 J ~ y-1 


(14) 


du p 

dt 


= — (u -u ) 

m s F 


(17) 


where R p and u p are position and velocity of particle 
along a streamline, u, is the velocity of the gas at the 
streamline location and D is a drag related term and is 
given by 


D= y 2 pPs c <ih- u p\ (18 > 

C d =-£L[l + J-Re\] for Re p < 1000 (19) 

Re p 6 y 

Q = 0. 424 for Re p > 1000 (20) 

where Re p is the Reynolds number based on the particle 
diameter. Equations (16) and (17) are integrated using a 
third-order variable step Runge-Kutta integrator. At 
each time step, u g is computed from equations (14) and 
(15). 


Initial Results and Discussion 


where a 0 is the stagnation speed of sound. The Mach 
number M is computed from the area Mach number 
relation 



t-J-M 2 ) ]<r+‘Mr-» ( 15) 


The area Mach number relation is solved by using the 
Newton's method. 


The model described above has been coded and used 
to predict the mass flux in a 5 meter radius 
hemispherical domain, resulting from the burst of a 
water jet from a n ominal Orbiter water dump. The 
diameter of the jet is 1.5 mm, the stagnation pressure is 
21.0 psia and the mass flow rate equals to 24.2 g/s. A 
vacuum water venting experiment was conducted at the 
AEDC 4- by 10-ft Research Vacuum Chamber (RVC) 
during 1983. Average particle sizes, particle number 
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densities and mass fraction were either measured or 
deduced from test data 4 * 5 . It was found that particle 
number density varies as 1/r 2 down the centerline of the 
plume. The number density in the polar direction has a 
bell shaped distribution from the centerline of the 
plume. These observations have been incorporated into 
the model. Two particle sizes woe chosen, large 
particle with mean diameter of 1.0 mm and small 
particle with mean diameter of 0.15 pm. The mass 
fraction of the large particle is 0.78; the mass fraction 
of the s m al l particle is 0.02; and the mass fraction of 
the gas is 0.20. Number density distribution of 
particles varies as 1/r 2 down the plume centerline; and a 
normal distribution function was used to describe 
number density distribution in the polar direction. The 
value of the number density at a reference point was 
scaled down from the AEDC test, such that the 
flow rate is equal to 24.2 g/s. The reference point is at 
1.0 cm down stream of the nozzle exit or the burst 
point. Standard deviation of the normal distribution 
function was set equal to either the bursting half angle 
or half of the bursting half angle. The former is 
designated as the la case and the latter as the 2a case. 
The results are then compared with the result obtained 
form an empirical model deduced form a video image of 
the STS-29 water dump 6 . Figure 6 shows the mass 
flux contours obtained from die empirical model. Figure 
7 shows the mass flux contours obtained from the la 
case. Figure 8 shows the comparison of mass flux at r 
= 5.0 m. This figure is plotted as mass flux versus 
polar angle, with 0 degree at the centerline of the 
plume. Figure 9 is an expanded view between 0 and 20 
degrees of Figure 8. Figure 10 shows the mass flux 
contours plot of the 2 or case. Figure 11 shows the 
comparison of mass flux at r = 5.0 m between the 2a 
case and the empirical model. Again, the figure is 
plotted as mass flux versus polar angle. Figure 12 is an 
expanded view between 0 and 20 degrees of figure 11. 
From these plots, we can see that fairly good 
comparison between the engineering model described in 
this document and the empirical model was obtained. 
The result from the empirical model falls somewhere in 
between the 7 a and the 2a cases. 

Model Improvements 

The current model for the Orbiter water dump is 
based on the semi-empirical results from the 
literature 2 * 4 * 5 . However, a mathematical model of the 
liquid flow into the vapor/ice cloud is desired. In 
general, as a liquid stream is discharged into a vacuum, 
due to the sudden drop of pressure, the stream becomes 
superheated and vapor bubbles may form inside the 
stream. As a bubble continues to grow, it eventually 
reaches a critical size which causes the jet to burst into 
water droplets and ice fragments. The fragments form a 
conically shaped plume. The burst distance and the 
characteristic cone angle are affected by the initial 


stream velocity, temperature and diameter. The 
technique of Smooth Particle Hydrodynamics (SPH) 
was chosen to build this model. The intent of the 
model improvement using SPH is to model the bubble 
formation and bursting 

Smoothed Particle Hy drodynamics 

SPH is a relatively new technique for hydrodynamic 
calcu l ations. It is a gridless Lagrangian method using a 
pseudo-particle interpolation method to compute 
smooth hydrodynamic variables, thus solving the 
Lagrangian equations of hydrodynamics. Each pseudo- 
particle has a mass, Lagrangian position and velocity, 
and internal energy, whereas other quantities are derived 
by interpolation or from constitutive relations. The 
major advantages of SPH are its relative ease of coding 
and its ability to handle irregular geometries and 
boundary conditions. A major disadvantage to SPH is 
it seldom "crashes" if something unphysical is done. 
Hence, one must continually ascertain that physical 
properties are conserved or realistic. Excellent reviews 
of the SPH method are given in Monaghan 10 * 11 and 
Benz 12 

A copy of an SPH code named SPHINX was 
obtained from the Los Alamos National Laboratory 
(LANL). This particular code had been tested and 
verified on a number of projectile impacts (see Figures 
13-16) and astrophysical problems 13 * 14 . However, the 
physical properties of radiative cooling, evaporative 
cooling, and surface tension were lacking in SPHINX. 
Work has been begun at JSC to implement these effects 
into an SPH formulation and add these models to 
SPHINX. Work has also been initiated at LANL to 
update SPHINX to make it more responsive to the 
Orbiter water dump problem. Unfortunately, the results 
of these efforts were not ready at the time of pre- 
conference publication. It is hoped that these results 
will be presented at the conference. 

Conclusions 

From the previous section, one can see that the 
result obtained from the engineering model described in 
this document agrees well with the empirical model 
deduced form the video image. The m^fn uncertainties of 
the engineering model are particle sizes, and number 
density distribution. To improve the model, further 
analysis of the available test data to derive correlation 
between particle sizes, number density, and distribution 
variance, versus initial stream condition will be 
conducted. Theoretical modeling of the bursting 
phenomenon will be accomplished to further improve 
the model. This update will focus on the use of the 
SPHINX code. 
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Bubble Growth in a Water Stream for Several Stream Temperatures 
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Figure 1 : Effect of Temperature on Bubble Growth 


Figure 3 : Composite Plot of Burst Distance versus 
Stream Temperature 


Bubble Growth an a Water Stream at Several Stream Velocities 
T = 29T It, D = 1.5 mm 
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Figure 2 : Effect of Stream Velocity on Bubble Growth Figure 4 : Burst Distance versus Stream Velocity 
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Figure 6 : Mass Flux Contours of Empirical Model 


Figure 8 : Comparison of Mass Flux at r=5m for the 
la case and the Empirical Model 
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Figure 1 1 : Comparison of Mass Flux at r=5m for the 
2c case and the Empirical Model 


Mass Flux at r = 5.0 m 



Figure 12 : Expanded View (0-20 degrees) of Figure 1 1 
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Figure 13 : SPHINX Simulation of the Impact of 

Aluminum Sphere on an Aluminum Plate 
t = 0.0944233 
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Figure 14 : SPHINX Simulation of the Impact of 

Aluminum Sphere on an Al uminum Plate 
t = 0.550503 ps 
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Figure 15 : SPHINX Simulation of the Impact of 

Aluminum Sphere on an Aluminum Plate 
t = 0.938149 ps 
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Figure 16 : SPHINX Simulation of the Impact of 

Aluminum Sphere on an Aluminum Plate 
t= 1.43971 ps 







